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In a recent Letter we determined analytically the right-handed quark mixing matrix in the min¬ 
imal Left-Right symmetric theory with generalized Parity. We derived its explicit form as a series 
expansion in a small parameter that measures the departure from hermiticity of quark mass matri¬ 
ces. Here we analyze carefully the convergence of the series by including higher order terms and by 
comparing with numerical results. We apply our findings to some phenomenological applications 
such as the production and decays of the right-handed gauge boson Wr, the neutrinoless double 
beta decay, the decays of the heavy scalar doublet, the strong CP parameter and the theoretical 
limits on the new mass scale from the K and B-meson physics. In particular, we demonstrate 
that the relevant coupling for the production of the Wr gauge boson at hadronic colliders and for 
the neutrinoless double beta decay equals its left-handed counterpart, within a percent. We also 
demonstrate that the stability of the theoretical lower limit on the Wr mass from the A'-meson 
physics is due to a partial cancellation of the external phases of the right-handed mixing matrix. 


I. INTRODUCTION 

The Left-Right symmetric theory [1] was put forward 
forty years ago as an attempt to understand the origin of 
parity violation in weak interactions through the spon¬ 
taneous symmetry breaking. The left-right (LR) sym¬ 
metry was chosen originally as a generalized parity ( V ), 
although it could also be a generalized charge conjuga¬ 
tion (C). Regardless of the nature of LR symmetry, the 
theory leads to a non-vanishing neutrino mass and to the 
seesaw mechanism [2-5] . The smallness of neutrino mass 
is related to the departure from the maximal violation of 
parity at low energies [2, 3]. Furthermore, the knowledge 
of neutrino masses and mixings allows to predict their 
Dirac Yukawa couplings and the associated physical phe¬ 
nomena [6]. 

Once the left-right symmetry is broken, the question is 
raised regarding the connection between the right-handed 
(RH) and left-handed (LH) quark and lepton mixing ma¬ 
trices. In the leptonic sector this connection goes away 
due to the seesaw mechanism which differentiate left and 
right handed neutrinos, while in the quark sector the sit¬ 
uation depends on the nature of LR symmetry. 

In the case of C , symmetric Yukawa couplings lead to 
symmetric quark mass matrices and thus to same left and 
right mixing angles, same magnitude but opposite signs 
of the Kobayashi-Maskawa (KM) phase and its right- 
handed analog, and to arbitrary external phases in the 
RH sector. The equality of mixing angles imply the same 
strength for production and hadronic decay rates of left 
and right-handed charged gauge bosons, while the new 
RH phases bring in new sources of CP violation. 

In the case of V, the situation is not that simple. In 
spite of having hermitian Yukawas couplings, the spon¬ 
taneous symmetry breaking of V introduces a complex 
phase which spoils the hermiticity of quark mass ma¬ 
trices, and a priori one would imagine no clear relation 
between LH and RH quark mixings. This soon became 
a burning issue, and the determination of the RH quark 


mixing matrix became over the years a kind of Holy Grail 
of the theory. Its explicit analytic form was only recently 
achieved [7] in the entire parameter space. In this longer 
sequel to our Letter we provide the missing technical de¬ 
tails and the higher order terms in the series expansion of 
Vr. We also discuss some important phenomenological 
applications and their consequences. 

A few words regarding the history of the attempts to 
determine the RH quark mixings. The first step in this 
direction was made in Ref. [8] through a numerical study 
in a portion of the parameter space. Some years later, 
Ref. [9] made an analytical study in the same approxima¬ 
tion and found the same hierarchical structure as in the 
CKM matrix. Finally, this result was established over 
the entire parameter space by combining analytical and 
numerical computations in Ref. [10]. What was missing 
until Ref. [7] was the exact equation for Vr , its analytical 
solution and a clear demonstration of the approximate 
equality of mixing angles of the left and right-handed 
sectors. 

A few years ago, a general low energy study of the RH 
quark mixing matrix was performed in Ref. [11]. More 
recently, it was argued that the near equality of the RH 
and LH mixing angles can be studied at the LHC with 
the b-quark tagging, and probed with a good precision 
at high luminosity [12], providing yet another motivation 
for our in-depth study of the RH mixing. 

The knowledge of the RH quark mixing is crucial for 
a number of phenomenological reasons, both at low en¬ 
ergies and at the hadronic colliders. 

First and foremost, the Majorana nature of left handed 
and heavy right handed neutrinos leads to neutrino¬ 
less double beta decay [13] through both left and right- 
handed gauge interactions [3, 14]. The RH contribution 
could easily dominate over the usual LH one, providing 
a natural example of new physics being responsible for 
this process [15]. 

Moreover, one can in principle observe lepton num¬ 
ber violation at hadronic colliders in the form of same 
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sign charged di-leptons, and probe directly the Majorana 
nature of heavy neutrinos through the so-called Keung- 
Senjanovic (KS) process [16]. Its potential discovery 
would allow for precise predictions [17] (see also [18]) 
for neutrinoless double beta decay and lepton flavor vi¬ 
olation; see e.g. [19] for a review. Detailed studies [20] 
support the feasibility of the KS process at the LHC and 
a possible probe of the LR scale all the way up to 5-6 
TeV; for a roadmap, see [21]. Recently the LHC has set 
the lower limit on the LR scale limit of roughly 3 TeV 
for a wide rage of RH neutrino masses [22]. It should be 
stressed that it is possible to determine the chirality of 
Wr gauge boson couplings, as discussed in the KS [23] 
and the hadronic channels [24]. 

While the hadronic colliders may be the best for a dis¬ 
covery of Wr due to their high energies, it is important to 
study also other collider possibilities. For a linear collider 
analysis we refer the reader to [25]. Recently, a study of 
the theory for the ep colliders was done in [26]. 

It was known for a long time that the small Kg — Kg 
mass difference leads to the lower limit on the LR scale 
in the minimal model [27], around 3 TeV (for the most 
recent study, see [28]). In the case of V, from the study 
of the electric dipole moment of the neutron, the limit 
was raised to about 7 TeV [29], without taking into ac¬ 
count the strong CP violation. Meanwhile new studies 
also improved the chiral perturbation results for this pro¬ 
cess [30, 31] and recently [31] claimed a limit (which de¬ 
pends though on the ultraviolet completion of the theory) 
of about 20 TeV by studying carefully the strong CP. 

It may be worth mentioning that the CMS reported 
recently [22] a 2.8 sigma excess in the KS process that 
can be interpreted as an indication of the LR symme¬ 
try [32]. It would require however the right-handed gauge 
coupling appreciably smaller that its left-handed counter¬ 
part, which would imply the breaking of LR symmetry at 
a very high scale. In what follows we work in the minimal 
model with single scale of symmetry breaking. 

The rest of the paper is organized as follows. In sec¬ 
tion II we discuss our main results, i.e., the derivation of 
the RH mixing matrix as a perturbation series in a small 
parameter which measure the hermiticity of the quark 
mass matrices. We show that mixing angles are well ap¬ 
proximated already at the first order; the differences be¬ 
tween RH and LH angles is in excellent agreement with 
the numerical results. The same is true for the KM phase 
and its RH counterpart, while the external phases in the 
RH mixing matrix need higher order corrections. We 
then discuss the issue of the convergence of the external 
phases. 

The phenomenological applications of our results are 
left for section III, where we discuss the production and 
decays of the RH charged gauge boson, the neutrinoless 
double beta decay, the decays of the heavy scalar doublet 
of the theory, the limits in the LR scale from the K and 
B meson physics and the computation of the strong CP 
parameter 9. 

Finally, in section IV we summarize our findings and 


offer an outlook for future work. 

In Appendix A we give the technical details necessary 
for the derivation of the analytic expression of the RH 
mixing matrix at the first order in perturbation theory. 
We discuss there how to deal with the square roots of 
matrices needed for this. In Appendix B we derive the 
the second and third order terms of the same expansion. 
In Appendix C, for illustrative purposes, we give exact 
values for the external phases in the case of vanishing 
mixing angles. 

Before we turn to the task of computing the RH mix¬ 
ing matrix, a comment is called regarding the question of 
potentially dangerous domain walls [33], the product of a 
spontaneous breaking of discrete symmetries. A natural 
way of avoiding them is of course inflation, but in a low 
scale theory it is not available. Another natural way [34] 
may be non-restoration of symmetries at high tempera¬ 
ture [35] , but even that may not work [36] . In any case, a 
tiny (even Planck scale suppressed) explicit breaking [37] 
suffices to get rid of domain walls. 


II. THE RIGHT HANDED QUARK MIXING 
MATRIX 


The LR symmetric theory studied here is based on the 
SU(2) l x SU(2)r x U(1)r-l gauge group augmented 
with generalized parity 'P: qg q r (for recent reviews, 
see [19]). Quarks and leptons come in LR symmetric 
representations 


Ql,r — 




(1) 


The Higgs sector consists of the following multiplets 
[3]: the bi-doublet $ and the SU(2)l,r triplets A l and 
A Rl where under generalized parity V\ <I> -f-> and 
A l tA A r . 

At the first stage of symmetry breaking, the neutral 
component of A^ develops a vev and breaks the original 
symmetry down to the SM one. The latter is in turn 
broken through the vevs of the neutral components of >I> 


(4>) = ?;diag(cos/3, — sin/3e “) (2) 


where v is real and positive and /3 < 7t/4, 0 < a < 2 tt. 

The quark Yukawa couplings in the minimal theory 
take the following form 


-L Y =q Tj (Y 1 3>-Y 2 <j 2 $>*<j 2 )q R + h.c. (3) 
which gives the following mass matrices for quarks 


( Mu \ =v ( ~ e “ S /3 UM 
\M d ) ^ e.~ la sp cp )\ Y 2 J 


(4) 


where cp = cos/3 and sp = sin/3. The underlying gener¬ 
alized parity V implies lrermitian Yukawa couplings 

yI 2 = y 1>2 


(5) 
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which in turn leads to the following relations between the 
up and down quark mass matrices 

M u - Ml = -is a t 2 p{e~ ia tpM u + M d ) (6) 

M d — M\ = is a t 2 p(M u + e ia tpM d ) (7) 

where s a = sin a, tp = tan/3, t 2 p = tan 2/3. The amount 
of the hermiticity of quark mass matrices is measured 
by the s a t 2 p parameter. From (7), by focusing on the 
third generation and ignoring its mixings with the first 
two, it is straightforward to obtain a rough upper limit 
s a t 2 p 2 mb/rnt- We justify it through an exact result 
found below. 

From the hierarchy of top and bottom quark masses, 
it is easy to see that M u is almost hermitian. Assuming 
small tan/3, Ref. [9] worked in the approximation of her¬ 
mitian M u and obtained a semi-analytical form of the RH 
quark mixing matrix. In what follows we derive an exact 


equation for Vr in the full parameter space and find an 
analytic solution as a series in the small parameter s a t 2 p. 
To set notation 


M u = U L m u Uh, M d = D L m d D^ R (8) 

where m q are diagonal matrices of positive quark masses. 
The left-handed CKM matrix Vl and its right-handed 
analog Vr are then given by 

V L = U}D L) Vr = U r Dr (9) 


as manifested in LH and RH charged gauge interactions 


C 


gauge — 


9 

V2 


{u L V L ty L d L + URVR^RdR) +h.c. 


( 10 ) 


We refer the reader to Appendix A for the details that 
go into the following leading order expression [7] 


( V R)ij = (V L )ij 


is a t 2 p 


tp{VL)ij + 


(VLm d Vl) ik (V L ) kj 
m Ui + m Uk 


{VL)ik(Vlm u V L ) kj 
m dk + m dj 


+ O{s 2 a t 2 20 ) 


( 11 ) 


This is not the only solution; the others are found 
through Vl —> S u VLS d and m qi —> s qi m qi , where S u = 
diag(s Ui ), S d = diag(sd i ) and s qi are ± signs. This is dis¬ 
cussed at the end of the Appendix A. Notice that even 
when quark mass matrices are hermitian, the LH and RH 
mixing matrices are not automatically equal due to the 
sign freedom of quark masses; instead one has the well 
known result Vr = S u VLS d . 

Equation (11) gives the RH mixing matrix Vr as an 
expansion in the small parameter s a t 2 p , which, for the 
sake of convergence, must satisfy 


\s a t 2 p\ < min 


m di + m dj 

(VlmMij 


9 ™± 
m t 


( 12 ) 


This agrees with what we estimated taking the third gen¬ 
eration and should not come as a surprise since the third 
generation CKM mixing angles are basically negligible in 
this context. 

When s a t 2 p is close to 2 mb/mt. one should include 
higher order terms, which can be found in the Ap¬ 
pendix B. This turns out necessary for the phases, while 
the values of the RH mixing angles obtained from (11) 
agree very well with the numerical solution of the exact 
equation (All). 

The first term in the bracket in (11) is quite small, 
while the third term generally dominates over the sec¬ 
ond due to large mass ratios of the top quark mass over 
the down quark masses (basically due to large mt/nib), 
and to some degree, the ratio of charm quark mass over 
m d and m s . It is noteworthy that the above results are 
parametrization independent. 


It is evident from (11) that a diagonal Vl implies a 
diagonal Vr at this order. We will see that this fact per¬ 
sists at higher orders and, moreover, it is clearly a solu¬ 
tion of the exact equation (All). The difference between 
RH and LH mixing angles is controlled by the smallness 
of the CKM angles and the smallness of the parameter 
s a t 2 p, as we discuss amply in what follows. 

A comment is in order. Our results, such as formula 
(11), are obtained by using generalized parity V and are 
valid at their face value only at the scale of LR symmetry 
breaking or above. They can be safely applied at high 
energy; however, in order to use them at low energies one 
needs to take into account their scale dependence. We 
are interested principally in a LR scale accessible to the 
LHC which makes these effects tiny and we will ignore 
them in the rest of this work. 


A. RH angles and KM phase. 

We compute next the difference between the left and 
right mixing angles and the KM phase and its RH coun¬ 
terpart. To do this we use the following parametrization 
for a general 3x3 unitary matrix 

V R = diag(e^\ e^ 2 , e^) V{0§, S R ) diag(e^, 1) (13) 

where V(9‘1,5l) = Vl is the standard form of the left- 
handed CKM matrix used by the Particle Data Group. 
Besides the RH analog Sr of the KM phase Sl (from now 
on both called KM phases), Vr contains five external 
phases which cannot be rotated away since we used all 
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the phase freedom in defining the usual CKM matrix in 
the left sector. 

A straightforward computation from (11) gives the 
leading terms for the differences between mixing angles 


0R - Ol 

n23 _ n23 

R L 

Or 3 ~ Of 


rrit 

— S a t 2/3 -S23S13S5 

m s 

m t m s 

— S a t2/3 -Sl2Sl3S<5 

m b m b 

m t m s 

—S a t2/3 -Si2S23 s <5 

m b m b 


and similarly for the KM phases 


Sr — Sl 


S a t 2p 


w c c 23 + m t s 23 

m s 


(14) 

(15) 

(16) 


(17) 


where, for simplicity, we defined Sij = sindf, c, 7 = 
cos and ss = sincfL. These are leading order terms, 
more complete expressions for the angles were given 
in [7], and the exact first order terms can readily be ob¬ 
tained from (11). 

It suffices to change the signs of quark masses accord¬ 
ingly to get all the other solutions. As shown below, the 
absolute values of the mixing angle differences are quite 
stable under these transformations, while the KM phase 
difference varies somewhat. 

Notice that the angle differences vanish in the limit of 
CKM phase Sl going to zero. This follows from the fact 
that in this limit, the first order terms in s a t 2/3 in (11) are 
purely imaginary and thus affect only the phases. From 
the above formulas, it is evident that the angle differ¬ 
ences are very small, suppressed by small CKM mixings. 
Moreover, in the case of 2-3 and especially 1-3 mixing an¬ 
gles there is an additional suppression of a small quark 
mass ratio m s /m b that compensates for the large mt/m b 
factor. As if there was a conspiracy in nature to keep the 
symmetry between LH and RH mixing angles in a world 
with broken parity. 

In Fig. 1 we plot in red lines these first order results, 
and with blue dots the numerical solutions of the exact 
equation. The first order is an excellent approximation 
and the agreement between the two is manifest through 
the whole physical range of s a t 2 p■ Notice that the phase 
difference Sr — Sl is multiplied with the factor sin#! 3 . 
The reason is that the phases Sl and Sr are always ac¬ 
companied with sin Of and sin Of (which are practically 
the same), respectively. 

The crucial point, as we noticed, is that the differences 
of LH and RH mixing angles are always proportional to 
another small LH mixings, which control their smallness. 
From Fig. 1 it seems that this holds true in higher or¬ 
ders of perturbation in s a f 2/3 and the proof comes from 
a discussion of a non-realistic two generation situation. 
Indeed, from (11), the difference of left and right mixing 
angles is zero, since Vl is real. In [7] we showed that the 
remarkable equality of Ol and Or is actually exact in the 
two-generation case, which then guarantees the smallness 
of the mixing angles differences at all orders. 



sin a tan 2 


FIG. 1. The differences between the right and left handed 
mixing angles and the KM phases. The first order terms are 
given by red lines, the blue dots denote numerical solutions of 
the exact equation. The agreement is manifest in the entire 
physical region Satzp < 0.54. 


B. RH external phases. 


First, we compute the external phases from (11) and 
(13) 


— U}3 + S a t2/3 


m c cl 3 +m t s% 3 m d c\ 2 +m. s s\ 2 
TO, 2 TO,, 


mt 

UJ 2 — W3 ~ S a f2/3x- 

2777-5 


(18) 

(19) 


f 2 m c cl 3 +m t sl 3 2 m c c| 3 +m t s| 3 \ 
UJ4 — W3 — S a t2@ I C 12 -bS 12 - ~ -) 

V m s 2 m d ) 


2 m c cl 3 + m t sl 3 

^5 — ^3 — S a t23 C 12 --- 

Zrric 


( 20 ) 

( 21 ) 


The expressions above are somewhat more precise that 
what we gave in the previous Letter version of this work. 
Unlike the expressions for the mixing angles and the KM 
phases, the external phases depend strongly on the sign 
transformations that connect different solutions. The 
above formulas should be taken as an example with all 
positive signs. It is straightforward to get more precise 
and complete expressions for all the cases. There is one 
subtlety to keep in mind: in some cases sign changes 
make the phases start from 7 r and not from zero, but 
that is easy to figure out. 

We plot these phases in Fig. 2. Again, the first order 
results are shown in red, and the numerical results in 
blue. Notice that in this case the results start diverging 
for larger values s 0 f 2 ,3 > 0.03, which simply implies the 
need for higher order terms in (11), as discussed below. 
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FIG. 2. Right handed external phases. In red lines we give 
the first order expressions and with blue dots the numerical 
solutions of the exact equation; they start diverging only for 
larger values s a t 2 p > 0.03. 


C. The impact of sign alternatives 

The reader must be worried about the different signs 
that plague (11), due to sign freedom of quark masses. It 
is easy to see in the first order terms that only the signs of 
mixing angle differences are really sensitive to the signs 
of quark masses. The case of the KM phase difference 
is more subtle, due to the contribution from both charm 
and top quarks. 

In Fig. 3 we plot the impact of the quark mass signs on 
the absolute values of mixing angle and KM phase dif¬ 
ferences computed at the first order in s a t 2 @- We include 
all the possible sign variations, which results in band like 
spreads of these quantities. The stability of mixing an¬ 
gle differences is remarkable; for the KM phases there is 
some dependence on the signs manifested in two smaller 
bands, very near to each other. When discussing these 
quantities phenomenologically one can safely ignore dif¬ 
ferent signs, at least at first order in s a t 2 p- 

As can be seen in (18)-(21), the external phases depend 
quite sensitively on the signs of quark masses, except for 
u >2 and W 3 . This is to be expected, for sign changes are 
equivalent to phase changes from zero to 7 r. One has to 
live with this aspect of the theory when discussing the 
CP violating processes and it must be taken into account 
when setting limits on the LR scale. 


D. Higher orders: external phases 

As one can see from Fig. 2, first order expressions for 
the external phases (18)-(21) are a good approximation 
only for small s a t 2 p- We now address the question of the 
convergence of the external phases as one approaches to 
the upper bound of s a t 2 @- 


FIG. 3. First order of the absolute value of right handed mix¬ 
ing angles and KM phases differences. All signs are included 
and the stability of the result is manifest. 


Using the standard parametrization for Ur, we find 
the higher orders of the phase W 3 (other phases can be 
treated in a similar way). For this, it is enough to focus 
on the 3-3 element of Vr 


(VR)33 = e“Wos^ 3 cosf^ 3 (22) 

To get w 3 , we first expand (Ur)33 up to third order 
using the expressions given in the Appendix B. Keeping 
only dominant terms up to third order, one obtains 


(Ur)33 ~ C23C13 


1 ■ . m t /, 


5 23 


+ ^ 13 ) 


1 

2 




m t V m s 2 

2 mt) mt, 23 


(23) 


We can see that, with the exception of the first, the odd 
terms in the expansion are negligible. Thus, summing 
the series, one has to a very good approximation 


TYlt / 9 2 \ 

Sa ^ 2/S 2TO6 1 S23 S 13 ' 

( 

From (22) and (24) one gets 

sino;3 ~ —s a t2p (1 + s 23 + s is) (25) 

The existence of the solution is then guaranteed by (12). 

To get the rest of the phases, it turns out to be suf¬ 
ficient to substitute 003 in(18)-(21) by the corrected one 
given above. This takes into account only higher order 
corrections to the dominant mt/mb term; corrections to 
the rest of the factors is straightforward but we have not 
included them due to their negligible effect. The final 


( Ur ) 33 - C23 ci 3 exp - z arcsin 
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sin a tan 2y8 


FIG. 4. External phases with sintJ 3 = —SatipmtIQmb) cor¬ 
responding to higher order expression (green line) and the 
numerical solution of the exact equation (blue dots). The 
agreement with the numerical solution is excellent and covers 
the whole parameter space. 


result is shown in Fig. 4. Notice that the phase W 3 con¬ 
trols the range of validity of the whole solution of Vr, 
and it stops at ~ ±7r/2 when s a t 2/3 reaches its upper 
bound (12). The agreement with the numerical results is 
now manifest in the whole parameter space. Once again, 
we confirm that the perturbation expansion in the small 
s a t 2/3 works well already at first order and the conver¬ 
gence of the external phases is under control. 

Last but not least, for illustrative purposes, the reader 
is referred to Appendix C for the exact solution of the 
RH phases in the non-realistic case of vanishing mixing 
angles. 


III. PHENOMENOLOGICAL APPLICATIONS 

A. RH charged gauge boson at colliders 

The physics of the Wr charged gauge boson needs the 
knowledge of Vr as much as one needs to know the CKM 
matrix in the case W. The proximity of LH and RH 
mixing angles tells us that the Wr production at the 
hadronic colliders proceeds basically in the same manner 
as the production of the W boson. This is useful since 
both CMS and ATLAS [22, 38] in most of their searches 
for a generic W' assume identical production rates as for 
the W boson. 

More precisely, the production strength of Wr is pro¬ 
portional to 


most a percent; thus to a great precision the production 
rate of Wr is equal to that of Wl ■ Similarly, the decay 
rates into top and bottom, or di-jets in general are to 
an excellent approximation given by assuming the same 
left and right couplings. Needless to say, it is possible 
to have LR symmetry broken at large scales; in such a 
case one has to run down our expressions to the rele¬ 
vant energies. The predictions will then be modified in a 
calculable manner. 

On the other hand, the leptonic decay rates of LH 
and RH gauge bosons are completely different in general, 
due to the seesaw mechanism behind neutrino mass. In 
this case Wr can decay into two jets and two same sign 
charged leptons [16], an appealing possibility of direct 
observation of lepton number violation. 


B. Neutrinoless double beta decay 


The LR symmetric model offers naturally a new contri¬ 
bution [3] to the neutrinoless double beta decay (Ov/3/3) 
through the exchange of Wr and RH neutrinos, in com¬ 
plete analogy with the usual contribution of Majorana 
light neutrinos. It goes through the RH gauge currents 
and RH neutrino mass, depending crucially on both lep¬ 
ton and quark RH mixing matrices. The total 0v[3f3 rate 
can then be described by the effective mass parameter 
[17] 


\M? +N \ = \\(M v ) e 


+ 


KV/Onl 4 ^ 

|(Rl)ii | 4 M^ R 


fc 2 \ 

M N/ e , 


(27) 

where fc is a measure of the neutrino virtuality, M u and 
Mjv are the mass matrices of light and heavy neutrinos, 
respectively. As we mention above, the lepton mixing is 
not predicted by the theory, and its determination would 
require the eventual direct observation of lepton number 
violation along the lines of the KS process. However, it is 
reassuring to know that the RH quark vertex (Vr) n, as 
shown in (26), can be predicted with arbitrary precision, 
in the same manner as in the case of the Wr production 
at hadron colliders. 

An additional contribution to Ov/3/3 mediated by the 
RH doubly charged scalar is disfavoured by lepton fla¬ 
vor violating constraints [17]. The rate can also proceed 
through the Dirac masses of neutrinos but, for RH neu¬ 
trinos visible at the LHC, it amounts to a sub-leading 
effect [ 6 ]. Nevertheless, these sub-leading contributions 
also depend on the quark mixing Vr and their precise 
evaluation is thus equally feasible. 


C. The theoretical limits from K and B physics 


|(Vr)ii| — C12 


rrit 

1 + S a t 2/3 -S12S23S13S5 

m. 


(26) 


It is easy to see that |(Vr)h | 2 differs from |(Vl)ii | 2 by at 


Another nice example of the relevance of Vr are the low 
energy limits on the LR scale from the kaon and B-meson 
physics. They depend crucially on the RH mixings, and 
we will illustrate it on a few examples. 
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K meson physics. Let us discuss the K — K mass 
difference. As is well known, in this case the dominant 
effect comes from the charm quark in the box diagram 
with the Wr and Wr charged gauge bosons. Only the 
real part of the amplitude contributes and it is propor¬ 
tional to the following form of the RH mixing matrix 


R e[A LR (K A')]cxRe[(V L ) 2 i(F L *) 22 (F R ) 2 i(^) 22 ] 

2 2 

— Sd s s c l 2 s 12 


X 




m c c‘c 


23 


+ rn t s 23 ) 


2 


(28) 


masses of heavy gauge boson. Its interactions can be 
deduced from the Yukawa in (3) 


T _I" M d + e ia S 2 pM u 

- L y =-Ql - h --- H 

V VC 2 0 

_I" Md . , * M u + e ia S2pMd . * 

+q L -*cr 2 /i- ia 2 H 

V VC2/3 


ur 


d R + h.c. 
(32) 


where 


where we have used the formula for Vr up to second or¬ 
der (for higher orders see Appendix B). Alternatively, ne¬ 
glecting the small difference between left and right mixing 
angles one gets cf 2 sf 2 cos(w 4 ~W 5 ). Using (20) and (21) 
one arrives again at (28). Due to different sign options 
in Vr, the cosine can be either —1 (and remain stable 
for any s a t 2 p because of the cancellation of the first and 
second term in the rest of the difference W 4 — W 5 ) or lie 
between 1 and ~ 0.7, in this case depending also on the 
value of s a t 2 ( g. 

Notice an interesting fact that helps the stability of 
the limit on the Wr mass. The difference 0 J 4 — 0 J 5 is sub¬ 
stantially smaller that the individual 0 J 4 and W 5 phases 
as seen from ( 20 )-( 21 ), and we can see clearly how the 
limit on the scale actually emerges. 

In a similar manner one can compute the imaginary 
part of the amplitude, which we leave as an exercise to 
the reader. 

B meson physics. In this case the relevant quantity 
to compute to the leading order, due to the top quark 
dominance, is (see e.g. [ 10 ]) 

r di OC (VS)33(Vk)3i (29) 

where di stands for strange (i = 2 ) or down (i = 1 ) quark. 
Neglecting the tiny difference between the left and right 
mixing angles, one gets the leading contribution 

r d cx c 23 ci 3 (si 2 s 23 - ci 2 c 23 si 3 e Il5i? ) e“ 4 (30) 

r s oc -c 23 ci 3 ci 2 s 23 e“ 5 (31) 

with Sr, C 04 and w 5 given approximately by (17), (20) 

and ( 21 ), respectively at first order (to account for higher 
orders corrections one should use also (25)). Again these 
analytical results facilitate numerical studies. 


h = C/301 + e ia S/3<j)2, H = —e ia sp(f>i + cpfc (33) 

h and H are the doublets with and without a vev, respec¬ 
tively, and (j> 1 and </> 2 are SU(2) doublets with Y = —1, 
forming the bi-doublet <f> = (<^ 1 , *cr 2 </> 2 ). The light dou¬ 
blet, h, is effectively the SM one, while H stands for the 
heavy doublet. Since H has to weigh more than 10 TeV, 
they can be taken with an excellent precision to be mass 
eigenstates. The Yukawa interaction of H in the physical 
basis becomes 


r o-o ^V L m d V^ + e “s 2 / 3 m^ 

-Lh — u Ul - Ur 


+ H V *d L 


VC 2/3 

o*~ 7 — Vlm u V R + e ia s 2 pm d 


VC2/3 


d R 


+ h - Tl ™K+£^ 2£1'h± UR 

VC 2 p 

TJ+ _ m u V R + e ia S 2 pV L m d 

— H^ul ■ - - -or + h.c 


(34) 


VC 2/ 3 


The above interactions depend crucially on Vr just like 
the gauge interactions of the heavy gauge bosons. This 
is relevant for the low energy processes such as Kr — Kr 
mass difference and for the production and identification 
of the heavy doublet H at future hadronic colliders. 

Regarding I\ — K mass difference. The neutral com¬ 
ponent of the heavy doublet adds coherently to the box 
diagram and will not affect the stability of the bound on 
the Wr mass. The tree level amplitude is proportional 
to 

Re[A H (K K)] a Re[(U L t m u y R ) 2 i(U^m„V L ) 2 i] (35) 


D. The physics of the heavy scalar doublet from 
the bi-doublet 

It has been known for a long time that the second dou¬ 
blet in the bi-doublet, due to the large amount of flavor 
violation it induces in neutral currents, has to have a 
large mass on the order 10-15 TeV or so (see for exam¬ 
ple [10], [11]). This is not a problem since it gets its 
mass [39] from the large triplet vev, responsible for the 


which, due to gauge invariance [40], has the same flavor 
dependence as in the usual Wr — Wr box diagram. Using 
as before higher orders of Vr one gets the same leading 
term as in (28). 

Similarly, the production and decays of the heavy dou¬ 
blet depend strongly on the quark mixings Vr and, for 
small values of /3, it may dominate the interaction. No¬ 
tice an important fact: H° decays principally in bb , a 
very clear signature. The other decay rates, such as say, 
tt depend on f3. 















E. The strong CP parameter 

As argued originally in [41], generalized parity guar¬ 
antees the vanishing of the QCD strong CP parameter 
9 before the symmetry breaking. This makes the effec¬ 
tive 9 calculable through the knowledge of V R , as done 
recently approximately in [31]. As a particularly nice il¬ 
lustration of the power of our results, in what follows we 
calculate exactly the first order of the series expansion of 
9 in terms of the small parameter s a f 2) g. 

In terms of the matrices that diagonalize the up and 
down quark mass matrices, one has 

9 = arg det U\U R D ] l D r (36) 

Now defining, as in Appendix A, U u = U\U R , U d = 
D\D r , and using V R Ud = U U V R one gets (noting also 
that arg det Vl = 0) 

9 = arg det Uy V R (37) 

Using the explicit form of U u as a matrix square root 
from (A9), it gives 

9 = arg det [V R + *s a f 2/3 ( tpe~ ia V R + m~ l V L m d )] (38) 

This is an tree-level expression, depending mainly on 
quark masses and mixing matrices. It could serve to 
compute 8 with an arbitrary precision. We can expand 
(38) in the small parameter s a f 2 ^ and obtain the leading 
term 1 

9 = s a t 2 0 1 Retr(^rn~ 1 V L m d Vl - (39) 

This expression and the exact one above, disagrees 
with the approximate form presented in [31]. The dom¬ 
inant term, though, is the same and is given by 8 ~ 
and could have been readily guessed by ig¬ 
noring the small third generation CKM mixing angles. It 
is obviously huge unless s a f 2/ 3 is vanishingly small. This 
forces a large bound on the LR scale as argued in [31]. 

How solid is this argument? In order to answer this, let 
us digress for a moment and discuss the question of strong 
CP in the SM. The fact that, experimentally, the strong 
CP parameter must be small 9 < 10 -10 has been coined 
the strong CP problem. But is this really a problem? In 
other words, is the perturbative contribution to 9 much 
bigger than the experimental upper bound? The answer 
is no, in fact it is much smaller as we now remind the 
reader. 

It is true that the SM does not predict the value of 8 
and it is a puzzle why it is so small, in view of the fact 
that the weak CP violation parameter is much bigger. 
It is also true that 9 appears to be divergent in pertur¬ 
bation theory, however, not before the sixth loop [42] 


1 Using the expansion arg det [1 + ieM] = eRe(trM) + 0(e 2 ). 


and for any reasonable value of the cut-off it is negligi¬ 
bly small. For example, for A cu toff — Mpi anc k one gets 
9 ~ 10 -19 [42], orders of magnitude below the experimen¬ 
tal limit. In other words, the question of the strong CP 
violation is a question of high energy physics, not a ques¬ 
tion of the SM itself. More precisely, the large scale that 
sets the value of 9 decouples effectively from low energy 
processes, in the same way the scale of grand unification 
decouples from low energies, except from fixing the weak 
mixing angle [43] (of course, there is the new possibility 
of inducing a proton decay, but this is consistent with 
the decoupling in usual interactions). 

A nice example of new high energy physics is the 
Peccei-Quinn (PQ) symmetry [44]. It fixes 9 no matter 
how large its scale is, and for practical purposes it decou¬ 
ples from low energy phenomena, the exception being a 
highly sensitive physics of the axion, just as proton de¬ 
cay in grand unification. The fact that the PQ symmetry 
does the job in the LR theory without affecting low en¬ 
ergy physics has been emphasised already in [31]. The 
message from this is simple: one should not worry about 
9 in theories at low, or relatively low, energies, such as 
the SM or the LR symmetric theory. This is a philosophy 
we will succumb to, and simply imagine a new physics 
such as PQ symmetry that takes care of the strong CP 
issue. In other words, we will allow s a t 2i g to take any 
value below its strict phenomenological bound in ( 12 ). 
Independently of how this is settled, and independently 
of the precise value of 9 one could obtain, however, a 
lower limit on the LR scale by considering the electric 
dipole moments of atoms, such as mercury. Although 
well measured, its theoretical value is plagued by large 
uncertainties, so setting a reliable limit will have to wait 
some more. 


IV. SUMMARY AND OUTLOOK 

In a recent Letter we were able to elucidate the long 
awaited form of the RH quark mixing matrix V R in the 
minimal LR symmetric model augmented with general¬ 
ized parity. We found exact equations, valid in the entire 
parameter space, which allow for its numerical determi¬ 
nation. We gave there an approximate form for V R as 
a leading term in the expansion of a small parameter 
which measures the departure from the hermiticity of 
the quark mass matrices. Moreover, we argued in favor 
of the proximity between left and right mixing angles, 
using the important fact that in the two generation case 
9 R = 9l- The small CKM mixing angles, together with 
a small ratio between strange and bottom quark masses, 
then guarantee practically equal mixing angles in the re¬ 
alistic three generation case. 

In this companion paper we provided a more detailed 
and complete discussion, paying special attention to the 
convergence of the employed expansion by calculating the 
second and third order of the series. In the case of mixing 
angles, we find an excellent agreement with the numerical 



9 


results already at the first order, while in the case of the 
phases in general one needs higher orders. 

We also discussed here a number of phenomenologi¬ 
cal applications that depend strongly on the form of the 
RH quark mixing. Our study shows the importance of 
the explicit analytical and numerical knowledge of Vr. 
In the case of Wr we demonstrate that the collider pro¬ 
duction and signature in di-jet decays is perfectly well 
described with the usual SM couplings of the W boson, 
except the opposite chirality. Similarly, we determined 
the strength of the right-handed gauge currents which 
control the rates of neutrinoless double beta decay. We 
also show how the decay rates of the heavy scalar doublet 
in the bi-doublet can be studied with great precision now 
that Vr is made transparent. As a nice example of the 
application of our results we give an exact and a leading 
term expression for the strong CP parameter 9. More 
important, we gave an explicit expression for the depen¬ 
dence of K l — Kr mass difference on Vr and showed that 
the stability of the resulting limit on the Wr mass is due 
to a partial cancellation of external phases. 

One last comment. The determination of Vr , as we 


showed here, requires only the minimal quark Yukawa 
sector and it does not depend, at least not at the tree 
level, on the situation in the leptonic sector or on the 
details of the parity breaking at the high scale. In 
other words, it does need the Majorana seesaw picture 
of neutrinos and could as well apply to the original LR 
model [1]. After all, neutrinos could be Dirac parti¬ 
cles [48] , no matter how appealing we find the seesaw pic¬ 
ture. All that matters here is that our findings hold true 
in any LR symmetric model with the single bi-doublet <F. 
While we prefer the seesaw picture which makes the LR 
theory complete and predictive when it comes to neutrino 
mass, it is up to experiment to decide. 
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Appendix A: Leading order 


In order to determine Vr one has to deal with a matrix equation involving square roots of matrices [7]. For simplicity 
and illustration we begin with the following square root 

\J m 2 + ieA (Al) 

where m is a diagonal matrix, A a matrix, and e an expansion parameter. The problem is to express the above square 
root in terms of the elements of m and A for small e. There are different ways of attacking this problem; we will work 
in 2 dimensions, where a root of matrix itself has a simple analytical form. The resulting expression turns out to be 
valid for any number of dimension. 

Using the formula for square root in the 2x2 case we have (taking all signs positive to ease the presentation) 


v t J- 2 \fd 

where d = det(m 2 + ieA) and t = tr(m 2 + ieA). Expanding first 


Vd, = nii m2 

1 


1 + i - tr(m 2 A) 


Vt + 2y/d m i + to 2 

and using them in (A2) one arrives at 


1 - 


ie trA + mini 2 tr(TO 2 A) 
2 (nil + m^) 2 


to 2 + ieA ) = rriiSij + ie -—-hO(e 2 ) 

/ a mi + mj 


(A2) 


(A3) 

(A4) 


(A5) 


It is easy to check that this expansion is valid regardless the number of dimensions. Still, here we outline the 
general steps needed to find the square root of a matrix M of arbitrary dimensions. This is useful also when solving 
the equation numerically. The first step is to decompose M into 

M = XJX- 1 (A6) 


where J is called the normal form of M. Once J and X are found (see [47] for the details), one finds the root as 


Vm = xyVjy~ 1 x~ 1 


( A 7 ) 
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Here, the matrix Y parametrizes a set of continuous solutions and it is only present in some very special cases. 

Coming back to the square root in (Al), now working in any number of dimensions, the above decomposition allows 
a straightforward expansion in terms of the small parameter. The result being once again formula (A5). 

In what follows we outline the steps needed to determine V R . It will be useful to introduce unitary matrices U u and 
U d which become diagonal sign matrices when the corresponding mass matrices are hermitian. In the notation of (8) 

Uu = U{U Rl U d = d[d r (A8) 

Then from (6) and (7) one finds 


ml + is a t 2f 3 (tpe ia ml + m u V L m d V^j (A9) 


m 2 d - is a t 2 /3 (tpe ia m 2 d + m d Vlm u V R ^ (A10) 

One has an additional relation which arises from the definition of the mixing matrices 


U u = 


m. u 


U d = — 
m. d 


VlUcL = U U V R 

Together with (A9) and (A10) it allows the determination of V R in terms of V R , m u , m d , a and (3. 
one can expand the square roots in U u and U d in (A9) and (A10) in powers of s a t 2 p 


tp/o \ , 0 


( Uu)ij — (S u ^)ij + is a t 2 j3 I ^ (^u)ij + ^ 


Rno 


■ m v 


tp. 


( V R ) ij 


( U d )ij — ( S d )ij is a t 2 p I ( S d )ij + „ - 

\ ^ “r TUd 


+ 0 {s a t 2 p) 


+ U{s 2 a t\p) 


(All) 

. Indeed using (A5) 

(A12) 

(A13) 


where S qi = diag(s g J, m qi = s qi m qi and s qi are ± signs. Using these expressions in (All) one can find V R as a power 
series in s a t 2 /s 


Vr = Vr ] - is a -t 2 pV£ ] + • • • (A14) 

The explicit form of and can be read-off from (11), with V^ = V R up to sign matrices, and V^ 1 ' 1 the term 
in brackets. There are 2^ 2n_1 ^ independent solutions for n generations, due to the square root nature of (A9) and 
(A10). The rest is found through V R —> S u V R S d and m qi —> s qi m qi , where S u = diag(s Ui ), S d = diag(sdj with, as 
before, s qi being ± signs. 


Appendix B: Higher order terms 

Let us turn again to the two-dimensional prototype case. From (A2) it is straightforward to get higher order terms 


— \f m 2 + iemA ] = Si. 


Aij 




Aj/c m^Aki 


m d + mj (mi + m].)(mk + m.j)(mi + mj) 

^■iki RRk\ A/ El fc 2 TTlk-2 ( 1 


(mi + m kl )(m fcl + m k2 )( m k2 + mj)(m l + m :j ) \m t + m k2 mj + m k 


+ 0(e 4 ) (Bl) 


This can be shown to holds in general for matrices of any dimension. 

We can now proceed to determination of the second and third order of Vr in terms of s a t 2 g. For the sake of space 
and clarity, we present only contributions of the type m u /m d due to the large top to bottom quark masses ratio, since 
these are the only ones that can compensate for the smallness of s a t 2 p. As we have seen already at the first order, 
all other terms are sub-leading at this level. For this purpose it is enough we take U u ~ S u and compute only higher 
orders of U d . 
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Using (Bl) one has up to the third order 

iylm u v R ) 


(Ud)ij — ( Sd)ij 'is a t2(3' 


13 + s 2 / 2 


_^ t C Vim u V R )i k (mdVim u V R ) k j 

m di + rh dj a 2/3 ( m di + rh dk ) (m dk + m dj ) (m di + m dj ) 

1 


i .^3 [Vj^Tn u VR^ik\[jW j dM[ j TTi u Vji)kik‘2(jriT'dM[ j Tiri u Vji)i~ 2 j l ± 

a 2/3 (rh di +rhd kl ){m dkl + m dk2 ) (m dk2 + m dj ) (rh di + m dj ) l m di +rh d 


(B2) 


m dj + m dki J 

where S d = diag(s c j i ), rh di = s di m di and s di are ± signs. Expressing also V R as a series of powers of s a t 2 p 

V R = vf - is a t 20 V£ ] - (s a t 20 ) 2 \ 4 2) + i{a a t 2p fV_i 3) + 0(s 4 a t%) (B3) 

and using it together with (B2) in (All) one gets, after equating each order of s a t 2 p, the following recurrence relations 

(B4) 

{Vlm u V^) kj 


(V^h = (S u V L S d )i 


(Y^h - 0 S U V L ) 


ik _ 


m dk + m dj 

(•</( 2 )\ ~ (Q V 1 (^L m u^R ^)kij 

\V R )ij — \^u * Ljiki 


(B5) 


(Vlm u vj i 0 ) )k 1 kArn d Vtm. u vj i 0) ) k2 j 

{m dki +m dk 2 )(m dk2 + m dj )(m dki +m dj ) 


(V { r\o - 0 SuV L )i 


\ ™d kl + m dj 

(Vjlm u Vft ^)kij (V]l m uV R ■ > )fc 1 fe 2 ^)k 2 jV{V^m u V R ^)k 1 k 2 ( m dVL rrl ^R ^)k 2 j 


(B 6 ) 


+ 


m dkl + m d] ( m dki + m dk2 ) ( m dk2 + m dj ) (m dfci + m d .) 

{Vlm u V R ) )k 1 k 2 {m d Vlm u V^)k 2 k 3 {rri d Vlm u V^ ) )k 3 i ( 1 __ 1 


(m dkl + m dk2 ) ( rh dk2 + m dk3 ) ( rh dk3 + rh dj ) (m dki + rh d] ) l m dki + rh dk3 rh d] + rh dk 


After solving these equations one arrives at the final result 

( TYl u VL ) j 


(vf k = {v L )ij 


(v^h * (Vl) 


ik ~ 


m dk 


m dk + m dj 

(™ dkl +m d .){m dki +m dk2 ){m dk2 

( V R } )ij - {VL)ikA V LmuV L ) klk 2 (Vlm u V L ) k 2 k 3 (Vlm u V L ) k3j 

m dk (m dk m d - m dk m d ) 


m. dj ) 


(B7) 

(B 8 ) 

(B9) 

(BIO) 

(Bll) 


(m dki +m dk2 )(m dk2 +m dk3 )(m dki +m dk3 )(m dki +m d .)(m dk2 +m dj )[m dk3 +m dj ) 


The other solutions can be found by changing, in the above expressions, Vl —> S u VLS d and m qi —► s qi m qi with s qi 
being ± signs. We see once again that a diagonal Vl implies a diagonal for V R , as was found at the first order. 
This confirms that small differences between RH and LH mixing angles are protected by small CKM angles. This 
fact is supported by the exact result of the equality of mixing angles in the two-generation case. We emphasize the 
parametrization independence of the above expressions, which should simplify the task of applying our results to 
physical processes; one can opt for a preferred parametrization. 

A few words regarding the convergence of the expansion. We have already given the necessary condition in (12) 
when discussing the first order terms, but we wish to reassure the reader that the higher order terms respect it. In 
the above equations there are potentially large terms, that however by inspection are seen to be at most of order one 


max 


max 


m dk 


~ 1 


m dkl + m. dj 

™ dkl (m dki m dk2 - m d . m dk3 ) 


{m dki + m d .)(m dki + m dk3 ){m dk2 + m dj ) 


~ 1 


(B12) 

(B13) 


and thus the convergence is guaranteed when (12) is satisfied. The convergence of the series was also checked by the 
comparison with the numerical results. 
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Appendix C: Zero mixing angles 


For the sake of illustration and simplicity, we go to the limit of vanishing LH mixing angles, which imply also 
diagonal Vr. This is useful for the case of small mixings, such as 6*23 and $13 and gives an idea of what happens in 
the general case with non-zero angles. From (All), one then finds Vr diagonal with 


±\/1 — S a t 2 @ 


( Vr)u = 


mz . - mj 


2 rrid m u 


- is a t20- 


2m d . m n 


1 is a t 2 p tpc 


The solution is valid when 


\s a t2p\ < 2 


I m di m Ui \ 

\ m li - m %\ 


(Cl) 


(C2) 


where one should take clearly the smallest value. Applied, for example, to the 2-3 generation case, this would give 
again roughly the limit 2mj/mj found before. In this case, by expanding in s a t 2 p one gets (ignoring corrections up 
to m b /m t ) 


{Vr)33 — ± 




1 

is a t 2i 3- 


m t 

m b 


(C3) 


Since in our convention (Vr) 33 = e*“ 3 , this leads to 


sin w 3 ~ -s a t 2 p 


1 m t 

2 m b 


As we show in (25), one gets basically the same result in the limit when the mixing angles are small. 


(C4) 
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